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Abstract 


Friction  of  graphite  and  diamond  surfaces  against  a  sharp  silicon  nitride  tip  was  measured  using  a 
friction  force  microscope  (FFM).  Atomic-scale  friction  of  a  freshly  cleaved  highly  oriented 
pyrolytic  graphite  (HOPG)  exhibited  the  same  periodicity  as  that  of  corresponding  topography. 
However,  the  peaks  in  friction  and  those  in  corresponding  topography  profiles  were  displaced 
relative  to  each  other.  Using  Fourier  expansion  of  the  interaction  potential,  we  have  calculated 
interatomic  forces  between  the  FFM  tip  and  graphite  surface.  We  have  shown  that  the  variations 
in  atomic-scale  friction  and  the  observed  displacement  between  the  peaks  in  friction  and  those  in 
corresponding  topography  can  be  explained  by  the  variations  in  interatomic  forces  in  the  normal 
and  lateral  directions.  At  large  scan  sizes  (50  nm  x  50  nm  or  larger),  the  variation  of  friction  for 
graphite  and  a  single-crystal  (Ila)  diamond  was  found  to  follow  the  local  slope  of  the  sample 
surface,  suggesting  that  a  ratchet  mechanism  is  operative  in  microscale  friction. 
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I.  INTRODUCTION 


Understanding  of  adhesion,  friction,  wear,  and  lubrication  of  interfaces  ("tribology")  is 
crucial  in  optimizing  the  performance  of  machinery  that  require  relative  motion  *'2.  Many 
physical  and  chemical  processes  affect  the  tribological  behavior  of  an  interface  in  relative  motion 
and  makes  the  understanding  of  the  underlying  mechanisms  difficult.  The  advent  of  new  tools- 
atomic  force  microscope  (AFM)  (first  developed  by  Binnig  et  al.^  in  1986)  and  its  modifications 
allow  the  study  of  contacts  on  a  micro-  to  atomic  scale.  By  modifying  an  AFM,  Mate  et  al.'* 
developed  a  friction  force  microscope  (FFM)  to  measure  atomic-scale  friction  of  a  shaip  tungsten 
tip  (with  a  radius  of  about  30  nm)  in  contact  with  a  test  surface,  which  is  a  convenient 
idealization  of  single  asperity  contact.  The  AFMs  and  FFMs  have  the  necessary  sensitivity  to 
allow  the  application  of  forces  weak  enough  (<  1  nN)  not  to  dislodge  atoms  from  their  sites 
during  contact.  Since  the  successful  measurements  of  friction  force  by  Mate  et  al.^,  several 
groups  have  combined  measurements  of  normal  and  friction  forces.  These  designs  have  used 
capacitance  measurements^,  laser  beam  deflection^-^  and  optical  interferometry®.  Now, 
commercial  instruments  are  available  which  allow  tiie  simultaneous  measurements  of  surface 
topography  and  normal  and  friction  forces. 

In  their  pioneering  study.  Mate  et  al.^  measured  atomic-scale  friction  of  a  tungsten  tip 
sliding  on  the  basal  plane  of  a  single  grain  of  highly  oriented  pyrolytic  graphite  (HOPG).  Later 
the  same  group  (Erlandsson  et  al.^)  measured  the  atomic-scale  friction  between  cleaved 
muscovite  mica  and  a  tungsten  tip.  For  both  graphite  and  mica,  these  authors  reported  that 
friction  between  a  tungsten  tip  and  the  test  sample  surfaces  exhibits  the  periodicity  of  the  sample 
surface  structure.  They  discussed  atomic  features  in  the  friction  force  in  terms  of  atomic-scale 
stick-slip.  The  magnitude  of  average  friction  force  varied  linearly  with  load  with  an  atomic-scale 
coefficient  of  friction  of  0.01  and  0.09  for  graphite  and  mica,  respectively.  Germann  et  al.*® 
measured  the  friction  of  a  diamond  tip  on  diamond  (111)  and  (100)  surfaces  in  ultrahigh  vacuum 
(UHV).  The  friction  force  for  (100)  surface  (with  a  step-like  structure)  exhibited  periodic 
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features.  They  suggested  that  these  features  were  associated  with  a  2x1  reconstructed  surface 
convoluted  over  an  asymmetric  tip  shape.  The  (111)  surface  showed  features  which  could  not  be 
simply  related  to  the  surface  structure.  Hirano  et  al.*  •  used  mica  instead  of  a  sharp  tip  and 
measured  the  friction  between  a  cylindrical  and  a  flat  single-crystal  muscovite  mica  samples. 

They  reported  that  frictional  forces  were  anisotropic  with  respect  to  the  lattice  misfit  angle 
between  the  two  mica  crystals.  To  date,  a  clear  understanding  of  friction  mechanism  on  a  micro 
or  nano  scale  does  not  exist. 

Friction  force  images  can  also  provide  additional  information  which  can  be  used  to 
distinguish  different  chemical  species  in  multicomponent  organic  films  as  shown  by  Ovemey  et 
al.*2.  The  same  group  (Meyer  et  al.‘3.  i4)  have  also  studied  the  lubrication  behavior  of 
Langmuir-Blodgett  films  using  friction  force  microscope  and  observed  a  10  fold  reduction  in 
friction  for  hydrocarbon  covered  surfaces  compared  to  that  of  bare  silicon  substrates.  (Also  see 
O'Shea  et  al.*^). 

In  this  paper,  we  present  data  of  simultaneous  measurements  of  surface  height  and 
friction  force  profiles  of  graphite  and  diamond  surfaces  against  a  silicon  nitride  tip  on  micro-  and 
atomic  scales.  The  diamond  surface  did  not  show  distinguishable  atomic  features,  therefore 
atomic  scale  data  are  only  reported  for  a  graphite  surface. 

n.  EXPERIMENTAL 

We  used  an  AFM/FFM  from  Digital  Instruments  (NanoScope  III)  for  this  study.  The 
AFM/FFM  used  here  (Fig.  1)  can  provide  simultaneous  measurements  of  friction  force  and  surface 
topography.  The  sample  is  mounted  on  a  piezoelectric  tube  (PZT)  scaimer  which  can  precisely 
"scan  the  sample  in  the  horizontal  (x-y)  plane  and  can  move  the  sample  in  the  vertical  (z)  direction. 
A  sharp  tip  at  the  free  end  of  a  Si3N4  beam  is  brought  in  contact  with  the  sample.  A  laser  beam 
from  a  laser  diode  is  focused  onto  the  back  of  the  cantilever  near  its  free  end.  The  cantilever  is 
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tilted  downward  at  about  10  deg.  with  respect  to  the  horizontal  plane.  The  beam  is  reflected  off  the 
end  of  the  cantilever  and  is  directed  through  a  mirror  onto  a  split  photodetector  with  four 
quadrants.  Two  quadrants  [top  (T)  and  bottom  (B)]  of  the  detector  are  used  for  measurement  of  the 
surface  topography.  As  the  sample  is  scanned  under  the  tip,  topographic  features  of  the  sample 
cause  the  tip  to  deflect  in  vertical  direction.  This  tip  deflection  will  change  the  direction  of  the 
reflected  laser  beam,  changing  the  intensity  difference  between  the  top  and  bottom  photodetector 
(AFM  signal).  A  feedback  circuit  is  used  to  modulate  the  voltage  applied  to  the  PZT  scanner  to 
adjust  the  height  of  the  PZT,  so  that  the  cantilever  vertical  deflection  or  the  normal  force  (given  by 
the  intensity'  '^i^^ference  between  the  top  and  bottom  detector)  will  remain  almost  constant  during 
scanning.  Thus  the  PZT  height  variation  is  a  direct  measure  of  surface  roughness  of  the  sample. 

Details  of  the  friction  measurement  techniques  have  been  described  elsewhere*^.  We  used 
the  method  in  which  the  other  two  (L  and  R)  quadrants  of  the  photodetector  (arranged  horizontally) 
are  used  for  the  measurement  of  friction  force  being  applied  to  the  tip.  The  sample  is  scanned  back 
and  forth  in  a  direction  perpendicular  to  the  longitudinal  axis  of  the  cantilever  beam.  Friction  force 
between  the  sample  and  the  tip  produces  a  twisting  of  the  cantilever.  The  laser  beam  is  reflected 
out  of  the  plane  defined  by  the  incident  beam  and  the  beam  reflected  from  an  untwisted  cantilever. 
This  produces  an  intensity  difference  of  laser  beam  received  between  the  left  (L)  and  right  (R) 
quadrants  of  the  photodetector.  This  intensity  difference  (FFM  signal)  is  directly  related  to  the 
degree  of  twisting,  hence  to  the  magnitude  of  friction  force.  One  problem  associated  with  this 
method  is  that  any  misalignment  between  the  laser  beam  and  the  photodetector  axis  would 
introduce  error  in  the  measurement.  By  adding  the  average  of  the  two  FFM  signals  obtained  by 
scanning  the  sample  in  two  opposite  directions  and  dividing  by  two,  and  then  subtracting  this 
component  from  either  profiles,  the  misalignment  effect  can  be  eliminated.  Voltages  corresponding 
to  normal  and  friction  forces  can  be  converted  to  force  units  by  following  the  force  calibration 
procedures  developed  by  Ruan  and  Bhushan'^.  By  making  measurements  at  various  normal  loads, 
average  value  of  coefficient  of  friction  is  obtained  which  then  can  be  used  to  convert  the  friction 
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profile  to  the  coefficient  of  friction  profile.  Thus  any  directionality  and  local  variation  of  friction 
can  be  easily  measured.  Surface  topography  data  can  be  measured  simultaneously  with  friction 
data  and  local  relationship  between  the  two  profiles  can  be  established.  During  AFM  /FFM 
measurements,  typical  scanning  speed  was  about  200  nm/s  in  the  sliding  direction  and  2  nm/s  in 
the  direction  perpendicular  to  the  sliding.  Microfiriction  measurements  were  made  over  200  nm  x 
200  nm  and  50  nm  x  50  nm  scan  areas,  whereas  atomic  scale  friction  were  measured  over  2  nm  x  2 
nm  and  4  nm  x  4  nm  scan  areas.  The  atomic  scale  images  were  magnified  and  only  a  1  nm  x  1  nm 
portion  will  be  shown  for  clarity. 

A  highly  oriented  pyrolytic  graphite  (HOPG)  [rms  roughness  (a)  =  0.3  nm  over  an  1  pm 
x  1  pm  cleaved  area]  and  a  polished  single-crystal  natural  diamond  (Ua)  (a  =  8.6  nm  over  an  1 
pm  X  1  pm  area)  were  used  in  this  study.  Graphite  was  freshly  cleaved  for  atomic  scale 
measurement.  After  atomic  scale  measurements,  the  graphite  surface  was  roughened  (a  =  8.8  nm 
over  an  1  pm  x  1  pm  area)  using  a  sand  paper  (2400  mesh,  3-pm  grit  SiC).  The  sample  was 
subsequently  ultrasonically  cleaned  to  remove  debris  from  the  surface.  The  roughened  surface 
was  used  for  microfriction  measurements. 

III.  RESULTS  AND  DISCUSSION 
A.  Micro-Scale  Friction 

A  gray-scale  plot  of  the  topography  and  corresponding  micro-scale  friction  profiles  of 
roughened  graphite  for  a  200  nm  x  200  nm  scan  size  is  shown  in  Fig.  2  (top  and  bottom, 
respectively).  Also  shown  in  this  figure  (middle)  is  the  surface  slope  along  the  tip  sliding 
direction  (derivative  of  the  surface  topography)  (Fig.  2b).  Light  color  corresponds  to  high  points 
in  topography,  slope,  or  friction  and  dark  color  corresponds  to  the  corresponding  low  values.  It  is 
clear  that  the  friction  profile  does  not  resemble  the  corresponding  surface  topography  profile. 
However,  the  friction  profile  and  the  slope  of  the  corresponding  topography  look  remarkably 
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similar,  which  is  true  down  to  a  50  nm  x  50  nm  scan  (Fig.  3).  Similar  trends  are  obtained  for 
polished  diamond,  Fig.  4. 


The  resemblance  between  the  slope  of  surface  roughness  profiles  and  the  corresponding 
friction  profiles  can  be  explained  by  the  ratchet  mechanism  initially  proposed  by  Makinson*^  and 
later  discussed  by  Tabor’^.  We  consider  a  small  tip  sliding  over  an  asperity  making  an  angle  6 
with  the  horizontal  plane.  Fig.  6.  The  normal  force  (normal  to  the  general  surface)  applied  by  the 
tip  to  the  sample  surface  W  is  constant.  Friction  force  F  on  the  sample  varies  as  a  function  of  the 
surface  roughness.  It  would  be  a  constant  poW  for  a  uniform  and  smooth  surface  in  the  presence 
of  "adhesive"  ftiction  mechanism.  The  force  components  along  (S)  and  perpendicular  to  (N)  the 
local  surface  of  the  sample  (ascending  part  of  the  asperity)  at  the  contact  point  are, 

N  =  W  COS0  +  F  sin0  (Al) 

S  =  F  COS0  -  W  sin0  (A2) 

Since,  when  sliding  occurs  S/N  =/io  (coefficient  of  friction  without  roughness  effect),  we  obtain 
the  value  of  coefficient  of  friction  including  the  roughness  (ratchet)  effects 

\L\-  F/W  =  (po  +  tan6)/(l-po  tan0),  (A3a) 

where  fij  is  the  local  coefficient  of  friction.  For  a  small  /to  tan0,  Eq.  (A3a)  can  be  rewritten  as 
p  1  ~  Po  +  tan0  (A3b) 

indicating  that  the  asperity  term  is  additive  to  the  local  coefficient  of  friction.  Similarly,  on  the 
descending  part  of  the  asperity  where  tan©  is  negative, 

P2~  Po  +  tan© 

=  Po  -  ltan©|  (A4). 

For  a  symmetrical  asperity,  the  average  coefficient  of  friction  the  AFM  tip  experienced  in 
traveling  across  the  whole  asperity  is 
Pave  =(Pi+P2)/2 

~Po(l+tan2e)  (A5) 
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if  Ho  tan0  is  small.  Therefore  a  symmetric  asperity  yields  a  contribution  to  coefficient  of  friction 
that  is  proportional  to  ^cHanW  averaged  over  the  surface. 

Other  mechanism  which  also  has  a  relationship  between  riiction  and  surface  shape  is 
ploughing  mechanism.  According  to  Bowden  and  Tabor*,  the  ploughing  component  of  friction 
with  tip  sliding  across  an  asperity  is  fip  ~  tan0.  Since  in  the  FFM  measurements,  we  notice  little 
damage  of  the  sample  surface,  the  contribution  by  ploughing  is  expected  to  be  negligible  and  the 
ratchet  mechanism  should  be  the  dominant  mechanism  for  the  local  variations  in  the  friction 
profiles.  With  a  tip  sliding  over  the  leading  (ascending)  edge  of  an  asperity,  the  slope  is  positive, 
it  ia  negative  during  sliding  over  the  trailing  (descending)  edge  of  the  asperity.  Thus,  friction  is 
high  at  the  leading  edge  of  asperities  and  low  at  the  trailing  edge.  The  ratchet  mechanism  thus 
explain  the  correlation  between  slopes  of  the  roughness  profiles  and  friction  profiles  observed  in 
Figs.  2  to  4.  We  note  that  in  the  ratchet  mechanism,  FFM  tip  is  assumed  to  be  small  as  compared 
to  the  size  of  asperities.  This  is  valid  since  the  typical  radius  of  curvature  of  the  tips  is  about  30 
nm.  The  radius  of  curvature  of  the  asperities  of  the  samples  measured  here  (the  asperities  that 
produce  most  of  the  friction  variation)  is  found  to  be  typically  about  100-200  nm  which  is  larger 
than  that  of  the  FFM  tip. 

B.  Atomic-Scale  Friction 

The  atomic-scale  topography  and  friction  profiles  of  freshly  cleaved  HOPG  are  shown  in 
Fig.  5  (a).  It  is  clearly  seen  that  friction  profile  exhibits  the  periodicity  of  the  graphite  surface 
structure.  However,  the  high  points  in  the  topography  and  the  high  friction  points  are  shifted 
relative  to  each  other.  This  can  be  seen  more  clearly  in  Fig.  5  (b).  Atomic  features  could  not  be 
resolved  for  diamond  and  the  corresponding  friction  data  for  diamond  are  not  presented  here. 

We  examine  the  atomic-scale  variations  of  friction.  Our  primary  purpose  is  to  understand 
the  spatial  shift  between  peaks  of  friction  and  those  of  the  corresponding  topography.  This  shift 
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reflects  the  spatial  displacement  of  maxima  of  normal  force  profile  relative  to  those  of  frictional 
force  profile  for  a  given  tip-sample  distance.  These  forces  can  be  calculated  from  the  derivative 
of  the  interaction  potential  (energy)  between  atoms  of  the  FFM  tip  and  those  of  the  sample 
surface.  We  will  thus  seek  an  analytical  expression  for  this  potential  followed  by  expressions  for 
the  normal  and  lateral  (friction)  forces.  Since  our  later  discussion  will  deal  with  the  topography 
of  graphite,  we  first  examine  one  of  the  relevant  features  of  topographic  images  of  graphite  - 
"asymmetric  atomic  structure". 

B.l.  Asymmetric  atomic  structure  and  topographic  images  of  Graphite 

We  first  show  in  Fig.  7  the  crystal  structure  of  (0001)  graphite  surface  which  consists  of 
multiple  six-membered  carbon  rings  ("benzene  rings").  Each  benzene  ring  consists  of  adjoining 
carbon  atoms  spaced  0. 1415  nm  apart  (C-C  distance).  Each  carbon  atom  has  three  neighbors 
arranged  at  the  apices  of  an  equilateral  triangle.  Three  of  the  four  valence  electrons  in  each 
carbon  atom  are  locked  up  with  the  tight  covalent  (sigma)  bonds  and  the  remaining  valence 
electron  is  delocalized  and  free  to  wander  over  the  surface.  Thus  adjacent  layers  have  only 
secondary  (van  der  Waals)  attraction  to  each  other.  Adjacent  layers  [(0001)  planes]  are  separated 
by  0.3354  nm  and  are  staggered  as  shown  in  Fig.  7.  Because  of  the  staggering,  only  one-half  of 
the  carbon  atoms  have  atoms  directly  above  and  below  (A  atoms)  in  the  adjacent  planes.  The 
other  half  (B  sites)  atoms  do  not  have  atoms  above  or  below;  they  fall  above  the  center  of  the 
hexagon  of  the  plane  below.  The  H-sites  (hollow  sites)  are  in  the  center  of  the  hexagon  (do  not 
have  carbon  atoms)  in  the  top  layer  and  have  atoms  directly  below  in  the  second  plane 
(corresponding  to  B-sites  in  the  second  plane).  The  7-atom  hexagon  connected  by  heavy  broken 
lines  in  Fig.  7b  (one  atom  is  at  the  center  of  the  hexagon)  represents  the  location  of  B  atoms  (to 
be  discussed  later). 

In  few  STM*^  (scanning  tuimeling  microscopy)  and  AFM  images^o,  every  atoms  in  the 
graphite  surface  have  been  observed.  However  in  most  STM21-24  AFM^^-Z?  measurements 
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of  graphite,  only  every  other  carbon  atoms  on  the  su.face  were  observed.  The  other  carbon 
atoms  were  either  not  as  prominent  or  "missing"  completely.  Thus  the  STM  and  AFM  images 
consist  of  benzene  rings  having  a  distance  between  two  nearest  prominent  atoms  of  about  0.246 
nm  with  an  additional  carbon  atom  at  the  center  of  the  ring  (Fig.  7b).  This  loss  of  trigonal 
symmetry  (between  the  A  and  B  types  of  atoms)  is  what  is  called  "asymmetric  atomic  structure" 
of  graphite  surface.  In  addition  to  the  "asymmetric  atomic  structure",  an  enormous  apparent 
height  ("corrugation  height")  of  atoms  above  the  centers  of  the  carbon  rings  (hollow  points  H, 
"holes")  has  also  been  reported.  For  example,  for  STM  images,  three  hills  and  three  saddle 
points  were  seen;  a  corrugation  between  hills  and  holes  as  high  as  2.4  nm  was  observed^^.  28. 

For  AFM  images,  a  corrugation  distance  of  0.02  nm  between  the  holes  and  the  highest  sites 
(designated  as  B  sites,  to  be  discussed  later)  and  0.002  nm  between  the  B  sites  and  the  second 
type  of  sites  (designated  as  A  sites)  was  observed^^. 

A  number  of  theoretical  papers  have  appeared  to  explain  the  STM*^-  23-24. 28-33  and 
AFM22.33-37  images  of  graphite  and  the  associated  phenomena.  Tomanek  et  al.^^*  3i  have 
calculated  the  electronic  structure  of  graphite  using  ab  initio  pseudo  potential  local  orbital 
method  and  concluded  that  the  asymmetry  is  purely  due  to  electronic  effect.  The  Fermi  surface  of 
graphite  collapses  to  a  single  point  in  the  surface  Brillouin  zone,  localized  at  A-  and  B-  sites. 

The  electronic  states  on  A-sites  are  decoupled  from  those  on  B-sites.  However,  a  coupling  exists 
between  A-site  atoms  in  the  adjacent  plane  and  this  coupling  causes  a  dispersive  band  centered  at 
Ep  for  electronic  states  located  at  A-sites.  Whereas  such  a  coupling  does  not  exist  between  the  B- 
site  atoms  on  neighboring  planes,  leaving  an  undispersed  state  Ep  at  B-sites.  The  STM,  which 
scans  a  narrow  energy  region  below  Ep  detects  all  B-site  states  and  only  a  very  small  fraction  of 
A-site  states.  Therefore  B-site  atoms  are  more  "visible"  to  STM  than  A-site  atoms. 

Batra  and  Ciraci33  have  performed  self-consistent  charge-density  calculations,  where 
surface  relaxation  was  considered  with  three  layers  of  graphite  representing  the  (0001)  surface. 


They  argued  that  since  A-type  atoms  have  neighbor  atoms  directly  below,  a  weak  interaction 
between  them  produces  a  weak  chemical  bond  leading  to  a  slight  accumulation  of  charge 
between  the  planes  of  graphite  under  the  A-sites.  This  accumulation  may  be  at  the  expense  of  a 
slight  depletion  from  the  vacuum  side  of  A-site  atoms.  The  B-sites  have  no  atoms  below  and 
hence  have  a  slightly  higher  charge  density  in  the  vacuum  region  coming  out  of  the  B-sites 
compared  to  the  A-sites.  The  charge  density  asymmetry  between  A-type  and  B-type  atoms  is 
even  greatly  amplified  at  Fermi  energy  Ep  and  the  STM  pattern  might  appear  at  times  to  show 
only  B  types  of  atoms  which  forms  a  centered  hexagon  of  side  0.246  nm  (Fig.  8).  The  H  sites 
have  the  lowest  charge  density.  Thus  the  highest  corrugations  should  be  between  B  and  H  sites. 
Batra  and  Ciraci^^  also  predicted  a  corrugation  values  of  about  0. 1  nm  based  on  their  charge 
density  calculation  at  Fermi  surface. 

For  understanding  of  AFM  images,  interatomic  forces  present  between  the  tip  and  the 
graphite  surface  need  to  be  included.  Batra  and  Ciraci^^  and  Abraham  and  Batra^  (also  see 
Gould  et  al.^"^)  calculated  the  interatomic  forces  using  a  self-consistent  field  method  allowing  for 
a  full  relaxation  of  the  electronic  charge  distribution33  or  surface  relaxations^.  They  found  that  if 
the  tip  is  in  close  proximity  with  the  sample  (2.7  a.u.,  comparable  to  the  first  nearest  neighbor 
distance  in  graphite),  force  at  the  B  sites  was  highest  followed  by  A  sites  then  by  H  sites. 

(Forces  at  A  and  B  sites  are  found  to  be  quite  similar  in  magnitude,  however.)  Batra  and  CiraciS^ 
argued  that  the  force  difference  between  A  and  B  sites  may  be  due  to  the  excess  charge  density  at 
the  B  sites  since  the  bands  at  Ep  are  centered  on  B  atoms.  At  higher  tip-sample  spacing,  the 
relative  distribution  of  forces  changes.  At  4.7  a.u.,  the  force  at  the  H  sites  is  higher  than  that  at  A 
and  B  sites.  Thus,  depending  upon  the  tip  proximity,  different  images  can  be  observed.  With  a 
tip  in  close  proximity  with  the  samples,  only  B  atoms  may  be  more  visible  as  seen  in  AFM 
images.  It  should  be  noted  that  Batra  and  Ciraci's^^  calculated  force  difference-betw^n  A  and  B 
sites  is  very  small,  not  significant  enough  for  a  firm  distinction  between  the  two  types  of  atoms. 
In  principle,  both  A  type  and  B  type  atoms  should  be  visible  to  the  AFM  tip33. 
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Mizes  et  al.*^  and  later  Gould  et  al have  offered  another  interpretation  for  the  STM  and 


AFM  images  of  graphite,  respecti .  ,;ly.  They  have  shown  that  the  asymmetry  and  a  variety 
experimental  images  cou!  J  ^  e  reproduced  by  calculation  if  it  can  be  assumed  that  there  are  two 
atoms  at  the  end  of  the  STM  or  AFM  tip.  Each  tip-atom  can  image  the  surface  and  a  relative 
phase  shift  between  the  two  images  "seen"  by  the  individual  tip-atom  would  produce  a  resultant 
(experimental)  image,  whose  pattern  can  be  different  from  the  honeycomb  structure  of  graphite 
lattice. 


While  these  explanations  for  the  graphite  STM  and  AFM  images  seem  to  be  convincing, 
they  probably  have  not  provided  the  complete  picture.  For  example,  as  reported  by  Land  et  al.^s, 
the  asymmetry  (appearance  of  every  other  carbon  atom)  was  observed  in  STM  images  even  for 
single-layer  graphite  deposited  on  Pt.  As  pointed  out  by  Land  et  al.38,  though  electronic 
interaction  between  the  graphite  and  the  Pt  substrate  is  strong,  this  does  not  allow  a  distinction 
between  two  unique  types  of  carbon  sites.  Therefore  the  observed  "asymmetric  atomic  structure" 
of  graphite  cannot  be  fully  explained  by  Batra  et  al.VO*^^  or  Tomanek  et  al.'s^'**  3i  theories. 

Mizes  et  al.'s’^  and  Gould  et  al.'s^'^  theories  do  not  have  the  difficulty  in  explaining  Land  et  al.’s^* 
data.  Yet  it  is  not  easy  to  verify  their  theories  since  the  precise  shape  of  the  STM  and  AFM  tips  is 
unknown.  As  we  shall  see  later,  there  may  be  other  factors  which  could  introduce  the 
"asymmetry  of  the  atomic  structure"  of  graphite. 

B.2.  Friction  force  images  of  graphite 

We  now  focus  on  the  spatial  shift  between  peaks  of  friction  and  those  of  the 
corresponding  topography  of  graphite.  As  mentioned  early,  this  shift  reflects  the  spatial 
displacement  of  maxima  of  normt  orce  relative  to  those  of  lateral  force,  and  these  forces  can  be 
calculated  from  the  spatial  derivative  of  the  interaction  potential  (energy)  between  atoms  of  the 
FFM  tip  and  those  of  the  sample  surface.  For  a  precise  analysis  of  the  interaction  force  between 
the  tip  and  sample  surfaces,  surface  relaxation  of  both  the  sample  and  the  tip  as  a  result  of  their 
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mutual  interaction  must  be  included.  An  exact  solution  to  this  complicated  interacting  system 
has  not  been  available  and  can  hardly  be  obtained  since,  first  of  all,  the  precise  shape  of  the  tip 
(physical  arrangement  of  the  outer  most  atoms  of  the  tip)  is  unknown.  In  most  theoretical 
calculations  the  tip  was  treated  as  "rigid"  or  as  a  single-atom  tip  with  no  relaxation,  only  the 
graphite  surface  was  allowed  to  relax2'^-28, 30. 33-37  With  this  and  some  additional 
approximations,  calculations  (such  as  the  one  using  continuum  elasticity  theory  with  parameters 
obtained  from  first  principles  total-energy  calculation  within  local-density  approximation35-37^ 
self-consistent  charge-density  calculation^®-  ^3)  have  provided  some  insight  towards 
understanding  the  observed  abnormalities  in  STM  and  AFM  images  of  graphite,  such  as  the 
observed  huge  corrugation  height23. 28  While  these  calculations  can  be  some  what  accurate 
depending  on  the  specific  problems  they  are  dealing  with,  they  are  performed  numerically  and 
are  time  consuming  particularly  if  tip  relaxation  is  allowed. 

Since  our  objective  here  is  to  understand  the  relative  spatial  shift  between  maxima  in 
topography  (or  normal  force)  and  those  in  friction,  a  simple  analytical  expression  for  the  normal 
and  lateral  forces  should  be  more  instructive  than  numerical  calculations.  These  forces  will  be 
calculated  from  the  spatial  derivative  of  the  interaction  potential  between  the  tip  and  the  sample. 
Approximations  are  necessary  in  deriving  such  an  expression.  The  simplest  approach  is  to  use 
Fourier  expansion  to  represent  the  interaction  potential.  Since  the  observed  topographical 
features  of  graphite  are  dominated  by  low  wave  vector  components,  it  is  sufficient  to  keep  only  a 
few  low  wave  vector  terms  in  the  Fourier  expansion  for  the  total  interaction  potential  which  leads 
to  the  desirable  simplicity.  We  will  make  two  assumptions;  (1)  the  interaction  potential  between 
the  tip  and  graphite  surface  (denoted  as  "total  interaction  potential")  equals  the  sum  of  the 
interaction  potential  between  individual  tip  atoms  and  carbon  atoms  of  graphite  (denoted  as 
"partial  potential"),  and  tlie  tip  atoms  interact  with  carbon  atoms  of  graphite  through  pair-wise 
potential  which  is  same  for  all  graphite  atoms;  (2)  the  pair-wise  interaction  potential  between  a 
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tip  atom  and  a  carbon  atom  can  be  approximated  by  Lennard-Jones  potential.  For  simplicity,  we 
will  consider  a  single  atom  to  represent  the  FFM  tip. 

The  periodic  structure  of  graphite  surface  (Fig.  9)  warrants  that  the  interaction  potential 
between  a  tip  atom  and  the  graphite  surface  is  also  periodic  with  the  same  periodicity  as  that  of 
the  graphite  surface  plane,  so  that  the  total  interaction  potential  u(r)  (r  is  the  position  of  the  tip 
with  respect  to  the  graphite  surface)  takes  the  form  of 

u(r+I)  =  u(r)  (Bl) 

with 

I=//a,+/2a2  (B2) 

where  /y  and  I2  are  integers  and  aj  and  a2  are  the  unit  lattice  vectors  in  the  surface  plane  (Fig.  9). 
The  periodic  nature  of  the  interaction  potential  [Eq.  (Bl)]  allows  u(r)  to  be  expressed  in  terms  of 
the  following  Fourier  expansion 

u(r)=  X  Wg(z)  exp(ig»t)  (B3) 

g 

where  X  is  the  two-dimensional  translation  vector  [in  the  (0001)  graphite  surface  plane],  z=lzl  and 
z  is  in  the  direction  perpendicular  to  the  graphite  surface,  g  is  the  multiples  of  the  reciprocal 
lattice  vectors  bj  and  b2: 

g=2jt(gybi  +  g2b2),  (B4) 

where  gj  and  g2  are  integers  and  bj  and  b2  are  defined  by 
ai.bi=l=a2.b2, 
ai.b2=0=a2.bt. 

Wg(z)  in  equation  (B3)  is  the  Fourier  coefficient  of  the  plane  wave  with  a  wave  vector  g.  There  is 
no  periodicity  in  the  direction  (z)  perpendicular  to  the  sample  surface  plane. 

Based  on  pair-wise  additivity^^  (assumption  1),  the  total  interaction  potential  u(r) 

between  the  tip  atom  and  the  graphite  crystal  can  be  regarded  as  the  sum  of  partial  interaction 
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potentials  (the  interaction  between  the  tip  atom  and  individual  carbon  atoms  of  graphite) 
u(r)=X  X  X  Uts(Za,X+l+mic)  (B5) 

a  /l,/2  k 

where  Uts(Za,X+l+ink)  is  the  partial  interaction  potential  between  the  tip  atom  and  the  carbon  atom 
in  the  ath  plane  in  the  graphite  crystal  located  at  a  vertical  distance  Za  from  the  tip  atom,  k 
indicates  the  kxh  atom  in  a  unit  cell  and  the  position  of  the  atom  in  the  plane  is  given  by  l+nik. 

Equations  (B3)  and  (B5)  are  different  expressions  of  the  same  physical  interaction.  Thus 
Wg(z)  in  Eq.(B3)  can  be  calculated  by  calculating  the  two  dimensional  Fourier  transformation  of 
Uts(Za.X+I+mk)  in  equation  (B5). 

w«(z)  =X  Wg(2a)  (B6) 

a 

with 

Wg(Za)  =-7  X  X  1  exp(-ig.x)  Uts(Za,X+l+mk)  dx 
^  lij2  k  A 

=— X  exp(ig.mk)[  exp(-ig.t)  UB(za,t)  dt  (B7) 

^  k  0 

where  A  is  the  smallest  area  representing  the  surface  periodic  structure^  (defined  in  this  paper  as 
"2-D  primitive  unit  cell"),  and  t=x+I+mk.  Equation  (B6)  means  that  the  coefficient  Wg(z)  of  the 
Fourier  component  with  wave  vector  g  is  the  sum  of  Wg(Za),  the  components  from  individual 
plane. 


Detailed  calculation  of  Wg(z)  requires  a  specific  function  for  Uts(r).  Steele^^  has  calculated 
Wg(Za)  using  the  Lennard-Jones  potential^  (assumption  2)  as  the  interaction  potential  between  a 


gas  atom  (tip  atom)  and  individual  atoms  in  a  solid  M,/r,), 


(B8) 
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where  r,  is  the  distance  between  the  tip  atom  and  the  ith  atom  in  the  solid.  The  parameter  which 
is  of  important  concern  here  is  the  value  a,  which  is  typically  between  0.2  to  0.4  nm^.  Using  this 
potential,  Steele^^  found  that 


Wo(z„)  =  e„^q( 


A 


^  4  ^ 


(B9) 


and 


w,(z„)=  [|l(-i-)=K,(g2„)-2o‘(^)'K,(gz„)J. 


'-a 


(BIO) 


where  q  is  the  total  number  of  atoms  per  2-D  primitive  unit  cell  (for  graphite,  q=2)  and  is  the 
modified  Bessel  function  of  the  second  kind.  Steele^^  also  mentioned  that  it  is  the  atoms  on  the 
surface  plane  which  make  appreciable  contribution  to  the  interaction.  This  is  because  both 
WoCza)  and  Wg(Za)  decrease  rapidly  with  the  increase  of  Za-  [Bassel  function  Kn(x)  goes  as 
(7t/2x)*'2exp(-x),  at  large  x.]  By  retaining  only  the  first  plane  (surface  plane)  in  the  summation  in 
Eq.(B6),  and  letting  Za  =  z,  we  have 

u(r)  ~wo(z)+XXLl  Wg(z)  exp[ig.(rak+i:)] 

g>0 


=C(z)+  SX|=i  D(g,  z)  exp[ig.(mk+x)] 

^>0 

27C 

where  g  is  the  magnitude  of  g,  g=lgl  =  — (gi^  +  g2^-gig2)*^  and 

a 


^  27c  .  2a 

C(z)=Wo(z„)  =  e^— q  (^^ 


4 

7  ^ 

a 


(BID 


(BI2) 


D(g,  z)=wg(za)=e^^  [£!l(^)*K,(gz„)-2a‘(^)^K,(gz„)].  (B13) 

There  are  two  carbon  atoms  per  2-D  primitive  unit  cell  in  graphite  basal  plane.  If  we 
choose  x=0  to  be  at  a  hollow  point  and  the  2-D  primitive  unit  cell  is  chosen  as  shown  in  Fig.  9, 
then  mic  for  the  two  atoms  in  the  2-D  unit  cell  indicated  in  Fig.  9  are  located  at  (1/3,  l/3)a,  (2/3, 
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2/3)a  (a  is  the  lattice  constant  of  the  2-D  unit  cell  tepresenting  the  distance  between  the  centers 

2jc 

of  nearest  hexagons).  Summation  over  nik  gives  a  value  of  2cos[— (g,  +  gj)].  Therefore 

u(r)  =C(z)  +2  X  D(g,z)exp(ig  •  x)cos[^(g,  +  gj )1  (B 14) 

«>0  ^ 


As  mentioned  early,  it  is  only  the  small  wave  vector  components  that  give  the  major 

contribution  to  the  interaction  potential  because  D(g,  z)  in  Eq.  (13)  decrease  rapidly  with  the 

increase  of  g.  To  the  first  order,  we  consider  only  the  first  six  wave  vectors,  for  which  Igl  =2n/a 

and  cos[2ic/3(gj+g2)]  =-1/2.  With  x  and  y  coordinates  chosen  as  shown  in  Fig.  9,  X  =x+y  can  be 

expressed  in  terms  of  a]  and  a2  with 
X 

X  =  -(ai-a2), 

a 


We  have 

u(r)  =  C(z)-D(27c/a,  z)  [2  cos  )  cos(— )  +  cos(-^)].  (B 15) 

aV3  a  a^l3 

Normal  and  lateral  forces  at  any  point  can  be  calculated  using  F=-V  u(r).  We  obtain 

27C  2itv  2‘Kx 

Fx  =  -D(27c/a,  z)  —  [2cos  )  sin  (— )];  (B 1 6) 

a  aV3  a 


Fy  =  -D(27c/a,  z)  [sin  (^)  cos  (— )  +  sin  (^)]; 

V3a  aV3  a  av3 


(B17) 


Fz  = 


dC(z)  dD(27C/a,z)  „  ^2^.  /4;iy^, 

^  [2  cos  (— ^ )  cos( - )  +  cos(— ^)]. 

dz  dz  a-43  a  a^3 


(B18) 


As  expected,  the  force  components  are  periodic  with  a  periodicity  a  in  the  x  direction  andVSa  in 

the  y  direction  as  shown  in  Fig.  9.  There  are  at  least  one  minimum  and  one  maximum  for  each 

force  component  within  each  2-D  primitive  unit  cell.  The  location  (xq,  yo)  of  the  extremes  of  each 

of  the  force  components  satisfies  the  conditions  that^* 

\d^Fi  d^Fi\ 


f  f 


and 


dxdy  dx^ 
d^Fi  d^Fi 


dy^  dxdy 
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<0, 


d^Fi  d^Fi 

where  i=x,  y,  or  z.  If  at  {xq,  yo)  -r-5-  <  0  and  — ^  <  0,  then  F,  has  a  maximum  at  this  point. 

dx  dy 

d^Fi  d^Fi 

Similarly,  F,  has  a  mimmum  if  —  >  0  and  —5-  >  0  at  (xq,  yo). 

ax  dy 


The  extremes  of  normal  and  lateral  forces  were  calculated  for  the  z  range  where 

<  0  and  D(27t/a,  z)>  0.  (In  the  region  where  ^  >  0  or  D(2jc/a,  z)<  0,  the 

dz  dz 

maxima  and  minima  of  corresponding  force  components  wiU  flop  their  positions.)  For  the 

normal  force  (F^,  the  minimum  is  at  (x=0,  y-0)  and  the  maxima  are  at  (x=0,  y=-7=  )  and  (x=0, 

-v3 

2a 

y=  -^)  within  the  2-D  primitive  unit  cell  indicated  in  Fig.  9.  The  positions  of  the  maxima  are 
precisely  where  the  carbon  atoms  are  located  at  and  the  position  of  minimum  is  at  a  "hollow 

■73 

point".  The  position  in  between  the  two  maxima  at  {x~0,  y=  — a)  is  also  a  special  point.  It  is 

2 

■73 

easy  to  see  that  if  the  x  coordinate  is  set  at  x=0,  then  Fj  has  a  minimum  at  y-—a  when  the  y 

■\/3 

coordinate  is  varied.  On  the  other  hand,  at  y  =  —a ,  F^  has  a  maximum  at  x=0  when  the  x 

73 

coordinate  is  varied.  Therefore  the  point  at  (x=0,  y=~a)  is  a  saddle  point.  Normal  force  Fj  at 

the  sad'^’'^  point  is  only  slightly  smaller  (by  10%)  than  that  at  the  two  maxima.  The  extremes  of 
the  lateral  forces  {Fx  and  Fj,)  were  calculated  and  are  summarized  in  Table  1  and  Fig.  9.  The 
extremes  of  Fx  along  the  x  axis  and  Fy  along  the  y  axis  are  also  shown  in  Fig.  9,  together  with 
detailed  variation  of  force  components  along  these  two  axes.  We  note  that  there  is  a  relative  shift 
between  maxima  of  lateral  force  and  those  of  normal  force.  The  shift  for  both  F*  and  Fy  is  in  the 
manner  that  the  lateral  force  is  to  push  the  FFM  tip  atom  away  from  the  points  where  the 
interaction  potential  is  the  largest. 


Examining  the  experimental  data  present^  in  Fig.  5  for  the  normal  force,  we  note  that 

only  one  maximum  was  observed  per  2-D  primitive  unit  cell,  i.e.,  only  every  other  atoms  were 

observed.  On  the  other  hand,  calculation  (within  the  first  order  approximation)  indicates  two 
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maxima  for  normal  force,  with  a  saddle  point  which  is  10%  smaller  than  the  maxima  in  between 
two  nearest  maxima.  (By  extending  the  calculation  to  second  order,  we  found  that  the  difference 
between  maxima  and  saddle  points  is  even  smaller.)  The  calculation  result  is  expected  since  we 
use  a  single  atom  to  represent  the  tip  and  only  the  lop  (surface)  layer  of  graphite  was  retained  in 
our  approximation.  By  including  more  layers  in  the  calculation,  modification  to  the  normal 
force  is  expected  to  be  small  since  the  force  decrease  rapidly  with  the  increase  of  the  distance 
between  the  tip  atom  and  the  underlying  layers.  This  is  consistent  with  Abraham  and  Batra's^'^ 
calculation  [where  six  (0001)  layers  were  considered].  They  found  that  the  (monoatomic)  tip 
experiences  very  similar  forces  when  it  is  placed  above  an  A  or  B  site.  Consequently,  by 
considering  only  the  normal  force  between  the  tip  and  the  graphite  surface,  a  monoatomic  tip 
would  detect  every  atom  of  the  surface  and  the  image  would  give  a  honeycomb  structure. 

The  discrepancy  between  the  calculated  and  the  experimental  (AFM/FFM)  data  for  the 
number  of  maxima  in  normal  force  profile  can  be  reconciled  by  assuming  that  there  are  two 
atoms  at  the  end  of  the  AFM  tip  using  Mizes  et  al.'s'^  and  Gould  et  al.'s^^  explanation  discussed 
previously,  or  assuming  that  there  is  a  graphite  flake  at  the  end  of  the  FFM  tip.  However,  the 
correctness  of  this  assumption  remains  to  be  verified  since  the  exact  shape  of  the  tip  is  unknown. 
We  note  that  in  all  reported  theoretical  calculations  of  topographic  imaging  of  graphite,  normal 
force  is  the  only  force  that  has  been  considered.  This  is  valid  if  the  tip-atoms  are  truly  "rigidly" 
fixed  to  the  tip  substrate.  However,  there  is  no  reason  to  believe  that  this  must  be  the  case. 
Intuitively,  if  a  sample  surface  (such  as  graphite)  is  allowed  to  deform  under  the  interaction  force 
between  the  tip  and  the  sample,  the  tip  surface  can  be  deformed  as  well.  In  the  case  of  a  single¬ 
atom  at  the  end  of  a  tip,  this  atom  may  relax  laterally  due  to  a  lateral  force  acting  from  the 
sample.  It  is  conceivable  that  the  atom  at  the  end  of  the  AFM/FFM  tip  can  actually  hop  around'^^ 
under  the  influence  of  lateral  (friction)  force.  For  example,  considering  an  FFM  tip  with  a  single 
atom  at  its  end  moving  across  a  graphite  surface  in  y  direction  as  shown  in  Fig.  9,  the  tip  atom 
experiences  a  lateral  force  towards  -y  direction  before  it  encounters  the  first  carbon  atom  of  the 
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indicated  2-D  primitive  unit  cell.  Lateral  force  is  zero  at  the  middle  of  these  two  carbon  atoms. 
After  the  tip  crosses  the  second  carbon  atom  of  this  2-D  unit  cell,  lateral  force  acting  on  the  tip- 
atom  is  strong  and  acts  towards  y  direction.  The  tip  atom  may  relax  from  one  position  to  another 
due  to  the  variation  of  lateral  force.  During  the  relaxation  the  normal  force  on  the  tip  will 
increase/decrease  slightly  which  will  be  compensated  by  the  sample  movement  (the  sample  will 
be  lowered/raised  to  maintain  the  normal  force  almost  constant).  Therefore  it  may  appear  that 
the  "real"  maximum  of  normal  force  as  seen  by  the  tip  atom  is  located  close  to  one  of  these  two 
carbon  atoms. 

Based  on  the  analysis  presented  in  this  section,  one  may  predict  the  relative  shift  between 
maxima  of  normal  force  and  those  of  friction  force  observed  in  Fig.  5.  K  we  assume  that  the 
FFM  tip  is  approximately  sliding  in  (x-0.16y  )  direction  (x  and  y  are  unit  vectors  along  x  and  y 
axes,  respectively),  and  that  the  normal  force  maxima  occur  at  the  location  of  the  second  carbon 
atom  in  the  indicated  2-D  primitive  unit  cell  and  its  equivalent  locations,  we  then  determine  a 
spatial  displacement  between  maxima  of  topography  and  those  of  friction  as  shown  in  Fig.  10. 
This  figure  looks  similar  to  Fig.  5b  from  the  experiment.  It  should  be  pointed  out  that  this  is  a 
simplified  analysis  involving  certain  assumptions.  More  detailed  calculation  involving  tip-atom 
relaxation  is  required  for  a  more  definitive  conclusion. 

Finally,  we  point  out  that  lateral  force  arising  from  atomic  interaction  is  a  conservative 
force  in  origin  which  would  yield  an  average  of  zero  friction  if  there  were  no  other  energy 
dissipation  mechanisms  involved.  A  complete  theoretical  treatment  relating  the  conservative 
lateral  force  to  true  frictional  force  has  not  been  available.  In  most  of  the  theoretical  calculations 
of  friction'^3'^,  it  was  simply  assumed  that  the  energy  that  is  required  for  a  tip  to  move  across  a 
potential  barrier  is  subsequently  (completely)  dissipated  into  heat.  Clearly,  additional  in  depth 
studies  are  required  to  understand  this  subject. 
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IV.  CONCLUSIONS 

The  mechanism  for  microscale  variation  of  friction  has  been  studied.  We  have  shown  that 
local  variations  in  friction  generally  follow  local  slope  of  sample  surface.  Ratchet  mechanism  is 
operative  in  this  scale.  We  have  also  developed  an  analysis  to  calculate  the  interatomic  forces 
between  atoms  of  the  AFM/FFM  tip  and  those  of  the  sample  surface  using  Fourier  expansion  of 
the  interatomic  potential.  We  have  found  that  atomic  scale  friction  is  related  to  these  forces. 
Maxima  of  normal  force  and  those  of  lateral  force  do  not  occur  at  the  same  location,  thus  there  is 
a  shift  between  the  peaks  in  the  friction  and  topography.  We  also  point  out  that  lateral  force  may 
play  an  important  role  in  topographical  imaging  of  graphite  surface. 
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Table  1  Location  of  extremes  of  force  components.  The  result  is  for  the  region  where 

<  0  and  D(2tc/a,  z)>  0  [a  is  the  lattice  constant  on  (0001)  graphite  plane 
dz 

and  is  also  the  distance  between  the  centers  of  nearest  hexagons  formed  of  carbon 
atoms  as  shown  in  Figs.  7  and  9]. 


Fx 

F. 

F. 

Maximum 

(l/4a,  ^a) 

(0. 1.474a) 

(0,  ^), 

2a 

Minimum 

(-l/4a,  ^a) 

(0. 0.258a) 

(0,0) 
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Figure  captions 


1 .  Schematic  of  atomic  force  microscope/friction  force  microscope  (AFM/FFM)  using  laser 
beam  deflection  method. 

2  Gray-scale  plot  of  surface  profile  (a=5.5  nm),  slope  of  the  surface  profile  (slope  profile) 
taken  along  the  sample  sliding  direction  ^mean  =  -0.09,  a=0.30),  and  friction  profile  (mean=l  nN 
0=2.5  nN)  of  a  200  nm  x  200  nm  scan  of  a  scratched  HOPG  for  a  normal  load  of  42  nN.  Light 
color  represents  high  points  in  topography,  slope  and  friction  profiles  and  dark  points  represent 
corresponding  low  values.  The  friction  and  the  slope  profiles  resemble  each  other. 

3.  Gray-scale  plot  of  surface  profile  (c=2.4  nm),  slope  profile  (mean  =  0,  a=0. 175)  and 
friction  profile  (mean=  InN,  0=1.7  nN)  of  a  50  nm  x  50  nm  scan  of  a  scratched  HOPG  for  a 
normal  load  of  42  nN,  showing  the  similarity  between  the  friction  and  slope  profiles. 

4.  Gray-scale  plot  of  surface  profile  (o=5.5  nm),  slope  profile  (mean  =  -0.006, 0=0.204) 
and  friction  profile  (mean=2.5  nN,  0=1.7  nN)  of  a  200  nm  x  200  nm  scan  of  a  polished  single¬ 
crystal  (Ha)  natural  diamond  for  a  normal  load  of  50  nN,  showing  the  similarity  between  the 
friction  and  slope  profiles. 

5.  (a)  Gray  scale  plot  of  surface  topography  and  friction  profile  of  a  1  nm  x  1  nm  area  of  a 
freshly  cleaved  HOPG  showing  the  atomic-scale  variation  of  topography  and  friction.  Higher 
points  are  shown  by  lighter  color,  (b)  Schematic  of  topography  and  friction  profiles  shown  in  (a). 
The  oblate  triangles  and  circles  correspond  to  maxima  of  topography  and  friction,  respectively. 
There  is  a  spatial  shift  between  the  two. 


27 


6.  Schematic  showing  the  effect  of  an  asperity  shape  on  local  variation  in  friction  in 
microscale  measurements. 

7.  (a)  Hexagonal  crystal  structure  of  the  (0001)  graphite  surface  showing  three  staggered 
layers,  (b)  Top  view  of  two  staggered  layers  (honeycomb  net)  shown  by  the  solid  (top  layer)  and 
dotted  (second  layer)  lines.  There  are  two  distinct  types  of  carbon  sites  in  graphite:  A-type  and 
B-type.  A-type  atoms  have  neighbor  atoms  directly  above  and  below  (solid  circles)  in  the 
adjacent  layers  and  B-type  atoms  (open  circles)  do  not  have  such  neighbors.  Points  "H"  on  the 
surface  are  "hollow"  points.  Based  on  analysis  presented  by  Batra  et  al^^.  and  Tomanek  et  al.^^, 
B-sites  atoms  are  more  prominent  than  A-sites  atoms,  thus  B-sites  atoms  are  generally  observed 
in  STM/AFM  imaging  (connected  by  heavy  broken  lines). 

8.  A  typical  STM  image  of  HOPG.  Only  every  other  carbon  atoms  in  the  surface  are 
"visible"  (bright  spots)  to  STM  tips,  the  other  half  atoms  are  "missing".  Similar  results  are 
obtained  for  AFM  (see  Fig.  5a). 

9.  Schematic  showing  the  orientation  of  the  x  and  y  axes  relative  to  the  surface  structure  of 
graphite.  Solid  lines  connecting  to  the  four  open  circles  form  one  of  several  possible  choices  of  a 
two-dimentional  primitive  unit  cell.  The  selected  location  of  2-D  unit  cell  simplifies  the 
calculation.  Carbon  atoms  are  located  at  the  vertex  of  hexagons.  Also  shown  are  locations  of 
extremes  of  normal  and  lateral  forces  between  a  graphite  surface  and  a  single  (FFM)  tip-atom  as 
well  as  force  variation  along  x  and  y  axes  based  on  a  theoretical  calculation. 

10.  Schematics  of  predicted  locations  of  maxima  in  normal  force  (topography)  (oblate 
triangle)  and  lateral  force  (circles),  assuming  that  due  to  lateral  force,  only  every  other  carbon 
atom  is  seen  in  the  normal  force  profile  (topography). 
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